Extrapolation Methods for the Dirac Inverter in Hybrid Monte Carlo

In Hybrid Monte Carlo(HMC) simulations for full QCD, the gauge fields evolve smoothly as a function of Molecular Dynamics (MD) time. Thus we investigate improved methods of estimating the trial solutions to the Dirac propagator as superpositions of the solutions in the recent past. So far our best extrapolation method reduces the number of Conjugate Gradient iterations per unit MD time by about a factor of 4. Further improvements should be forthcoming as we further exploit the information of past trajectories.Comment: latex file with espcrc2 styl

QCDF90: Lattice QCD with Fortran 90

We have used Fortran 90 to implement lattice QCD. We have designed a set of machine independent modules that define fields (gauge, fermions, scalars, etc...) and overloaded operators for all possible operations between fields, matrices and numbers. With these modules it is very simple to write high-level efficient programs for QCD simulations. To increase performances our modules also implements assignments that do not require temporaries, and a machine independent precision definition. We have also created a useful compression procedure for storing the lattice configurations, and a parallel implementation of the random generators. We have widely tested our program and modules on several parallel and single processor supercomputers obtaining excellent performances.Comment: LaTeX file, 8 pages, no figures. More information available at: http://hep.bu.edu/~leviar/qcdf90.htm

Towards a lattice calculation of the coefficients of the QCD chiral Lagrangian

We discuss a general strategy to compute the coefficients of QCD chiral Lagrangian by using the lattice regularization of QCD with Wilson fermions. This procedure requires the introduction of an effective Lagrangian for lattice QCD as an intermediate step in the calculation. The continuum QCD chiral Lagrangian can be then obtained by expanding the lattice effective Lagrangian in increasing powers of the external momenta. A suitable renormalization procedure is required to account for the chiral symmetry breaking introduced by the Wilson term in the lattice action. In anticipation of a numerical simulation, the lattice effective Lagrangian is computed analytically and investigated in the strong coupling and large $N$ limitComment: Talk presented at LATTICE96(chirality in qcd) , 3 pages, no figures. Latex file with espcrc2 styl

Monte Carlo Simulations of the SU(2) Vacuum Structure

Lattice Monte Carlo simulations are performed for the SU(2) Yang Mills gauge theory in the presence of an Abelian background with external sources to obtain information on the effective potential. The goal is to investigate the lowest Landau mode that, in the continuum one-loop effective potential, is the crucial mode for instability. It is shown that also in the lattice formulation this lowest Landau mode plays a very peculiar role, and it is important for the understanding of the vacuum properties.Comment: 3 pages, to appear in the Proceedings of Lattice 93, preprint BU-HEP-93-2

Chronological Inversion Method for the Dirac Matrix in Hybrid Monte Carlo

In Hybrid Monte Carlo simulations for full QCD, the gauge fields evolve smoothly as a function of Molecular Dynamics time. Here we investigate improved methods of estimating the trial or starting solutions for the Dirac matrix inversion as superpositions of a chronological sequence of solutions in the recent past. By taking as the trial solution the vector which minimizes the residual in the linear space spanned by the past solutions, the number of conjugate gradient iterations per unit MD time is decreased by at least a factor of 2. Extensions of this basic approach to precondition the conjugate gradient iterations are also discussed.Comment: 35 pages, 18 EPS figures A new "preconditioning" method, derived from the Chronological Inversion, is described. Some new figures are appended. Some reorganization of the material has taken plac

Polarizing Bubble Collisions

We predict the polarization of cosmic microwave background (CMB) photons that results from a cosmic bubble collision. The polarization is purely E-mode, symmetric around the axis pointing towards the collision bubble, and has several salient features in its radial dependence that can help distinguish it from a more conventional explanation for unusually cold or hot features in the CMB sky. The anomalous "cold spot" detected by the Wilkinson Microwave Anisotropy Probe (WMAP) satellite is a candidate for a feature produced by such a collision, and the Planck satellite and other proposed surveys will measure the polarization on it in the near future. The detection of such a collision would provide compelling evidence for the string theory landscape.Comment: Published version. 15 pages, 8 figure

Nonlinear Realization of Chiral Symmetry on the Lattice

We formulate lattice theories in which chiral symmetry is realized nonlinearly on the fermion fields. In this framework the fermion mass term does not break chiral symmetry. This property allows us to use the Wilson term to remove the doubler fermions while maintaining exact chiral symmetry on the lattice. Our lattice formulation enables us to address non-perturbative questions in effective field theories of baryons interacting with pions and in models involving constituent quarks interacting with pions and gluons. We show that a system containing a non-zero density of static baryons interacting with pions can be studied on the lattice without encountering complex action problems. In our formulation one can also decide non-perturbatively if the chiral quark model of Georgi and Manohar provides an appropriate low-energy description of QCD. If so, one could understand why the non-relativistic quark model works.Comment: 34 pages, 2 figures, revised version to be published in J. High Energy Phys. (changes in the 1st paragraph, additional descriptions on the nature of the coordinate singularities in Sec.2, references added

Chiral perturbation theory at O(a^2) for lattice QCD

We construct the chiral effective Lagrangian for two lattice theories: one with Wilson fermions and the other with Wilson sea fermions and Ginsparg-Wilson valence fermions. For each of these theories we construct the Symanzik action through order $a^2$. The chiral Lagrangian is then derived, including terms of order $a^2$, which have not been calculated before. We find that there are only few new terms at this order. Corrections to existing coefficients in the continuum chiral Lagrangian are proportional to $a^2$, and appear in the Lagrangian at order $a^2 p^2$ or higher. Similarly, O(4) symmetry breaking terms enter the Symanzik action at order $a^2$, but contribute to the chiral Lagrangian at order $a^2 p^4$ or higher. We calculate the light meson masses in chiral perturbation theory for both lattice theories. At next-to-leading order, we find that there are no order $a^2$ corrections to the valence-valence meson mass in the mixed theory due to the enhanced chiral symmetry of the valence sector.Comment: 25 pages, LaTeX2e; references adde

Effective Lagrangian for strongly coupled domain wall fermions

We derive the effective Lagrangian for mesons in lattice gauge theory with domain-wall fermions in the strong-coupling and large-N_c limits. We use the formalism of supergroups to deal with the Pauli-Villars fields, needed to regulate the contributions of the heavy fermions. We calculate the spectrum of pseudo-Goldstone bosons and show that domain wall fermions are doubled and massive in this regime. Since we take the extent and lattice spacing of the fifth dimension to infinity and zero respectively, our conclusions apply also to overlap fermions.Comment: 26 pp. RevTeX and 3 figures; corrected error in symmetry breaking scheme and added comments to discussio

Metallo-dielectric diamond and zinc-blende photonic crystals

It is shown that small inclusions of a low absorbing metal can have a dramatic effect on the photonic band structure. In the case of diamond and zinc-blende photonic crystals, several complete photonic band gaps (CPBG's) can open in the spectrum, between the 2nd-3rd, 5th-6th, and 8th-9th bands. Unlike in the purely dielectric case, in the presence of small inclusions of a low absorbing metal the largest CPBG for a moderate dielectric constant (epsilon<=10) turns out to be the 2nd-3rd CPBG. The 2nd-3rd CPBG is the most important CPBG, because it is the most stable against disorder. For a diamond and zinc-blende structure of nonoverlapping dielectric and metallo-dielectric spheres, a CPBG begins to decrease with an increasing dielectric contrast roughly at the point where another CPBG starts to open--a kind of gap competition. A CPBG can even shrink to zero when the dielectric contrast increases further. Metal inclusions have the biggest effect for the dielectric constant 2<=epsilon<=12, which is a typical dielectric constant at near infrared and in the visible for many materials, including semiconductors and polymers. It is shown that one can create a sizeable and robust 2nd-3rd CPBG at near infrared and visible wavelengths even for a photonic crystal which is composed of more than 97% low refractive index materials (n<=1.45, i.e., that of silica glass or a polymer). These findings open the door for any semiconductor and polymer material to be used as genuine building blocks for the creation of photonic crystals with a CPBG and significantly increase the possibilities for experimentalists to realize a sizeable and robust CPBG in the near infrared and in the visible. One possibility is a construction method using optical tweezers, which is analyzed here.Comment: 25 pp, 23 figs, RevTex, to appear in Phys Rev B. For more information look at http://www.amolf.nl/research/photonic_materials_theory/moroz/moroz.htm